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Abstract 



The inclusive production of the meson resonances ^"(770), /o(980) and /2(1270) in 
neutrino-nucleus charged current interactions has been studied with the NOMAD detector 
exposed to the wide band neutrino beam generated by 450 GeV protons at the CERN SPS. 
For the first time the /o(980) meson is observed in neutrino interactions. The statistical 
significance of its observation is 6 standard deviations. The presence of /2(1270) in neu- 
trino interactions is reliably established. The average multiplicity of these three resonances 
is measured as a function of several kinematic variables. The experimental results are com- 
pared to the multiplicities obtained from a simulation based on the Lund model. In addition, 
the average multiphcity of ^"(770) in antineutrino - nucleus interactions is measured. 

Key words: inclusive, resonance production, neutrino interactions 
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1 Introduction 



Hadron production in high energy neutrino-nucleon charged current (CC) inter- 
actions can be represented as a two-step process. The first step is the interaction of 
a W boson with a quark of the target nucleon changing it into a quark of a different 
flavour. The second step is the hadronization of the quarks during which coloured 
partons fragment into colourless hadrons. Perturbative QCD not being applicable 
to the soft process of hadronization, phenomenological models must be used. The 
most successful models are the string [1] and cluster fragmentation [2] models. 

A significant fraction of the particles in the hadron jets directly observed are 
produced through the production and decay of a large variety of very short lived 
particles, the so-called resonances. Consequently, the study of inclusive meson and 
baryon resonances is believed to reveal more directly the primary interaction mech- 
anism than the studies of stable particles (such as pions and kaons). Therefore, res- 
onance production provides a better ground for comparisons between models and 
experimental results. Inclusive resonance production has been extensively studied 
in e+e~ annihilation [3-8], photoproduction [9], electroproduction by muons [10], 
in hadronic interactions [11]. 

The role of the orbitally excited mesons, for example /o(980)(J^'^ = 0"'""'") 
and /2(1270)(J^^ = 2++), is of special interest in view of the possibly different 
production mechanisms. Both mesons have an orbital angular momentum equal to 
1 and a sum of spins of the valence quarks also equal to 1 resulting in a total spin 
of for the /o(980) and of 2 for the /2(1270). 

A number of nonstandard interpretations of the /o(980) state have been sug- 
gested. One approach regards the /o(980) as a four quark bound state [12], another 
as a two-kaon molecule [13]. Gribov [14] has proposed a new theory of QCD con- 
finement, in which the /o(980) plays the special role of a novel vacuum scalar state. 
Close et al. [15] noticed that the signature of Gribov's vacuum scalar state would 
be an enhanced yield in low multiplicity events. Robson [16] proposed an interpre- 
tation of the /o(980) as a scalar glueball. A recent OPAL study of the /o(980) at 
LEP [4] does not confirm the Gribov interpretation. All measured characteristics of 
/o(980) production in Z° decays are consistent with the /o(980) being a conven- 
tional scalar meson. However, an analysis of the WA102 data by Close and Kirk 
[17] selects the 'enigmatic' /o(980) as a glueball candidate together with /o(1500), 
/j(1710)and/2(1900). 

Bubble chambers used to be the only detectors where inclusive production of 
p and other resonances in neutrino interactions could be studied. In such detectors 
detailed information about the final state particles in each event can be obtained. 
The neutrino experiments [18-30] studied inclusive production of p°(770) and, in 
some cases, of /2( 1270) resonances. BEBC WA59 [20] studied p(770) (all three 
charge states), 7^(550), a;(783) and /2(1270). The resonance production properties 
were measured as a function of the mass of the hadronic system, of the charged 
hadron multiplicity and of various other kinematical variables. 

The main problem common to these studies was the small statistics (the total 
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of all experiments amounts to about 60000 v^CC events). The level of accuracy 
corresponding to this statistics did not necessitate taking into account the reflec- 
tions of other resonances (for the definition of reflections see Section 5). They were 
not accounted for in most of the previous neutrino experiments. Due to low statis- 
tics there were no attempts to look for the /o(980) meson in neutrino data. As to 
the other orbitally excited meson, BEBC WA21 [19] using 17750 v^p interactions 
found 615±226 /2(1270) mesons. But later, in the combined WA21 + WA25 + 
WA59 + E180 [31] data sample (58866 v^^CC events), no observation of /2(1270) 
was reported. 

All previous experiments reported an overestimation of the average multiplici- 
ties of p°(770) by the Lund model by 50-70%, while the shapes of the Lund distri- 
butions were in agreement with the experimental distributions. 

The NOMAD experiment in 4 years of data taking (1995-1998) has collected 
about 1.3 X 10^' v^CC interactions. In addition, a short antineutrino run has pro- 
vided about 32000 z/^CC events. The large statistics of the NOMAD data and the 
good quality of the event reconstruction make possible the study of processes that 
could not be studied in previous neutrino experiments. 

The article is organized as follows: section 2 describes the experimental setup 
including the NOMAD detector, the neutrino beam and the NOMAD simulation 
software; in section 3 the kinematical variables used in this analysis are defined; 
section 4 and section 5 describe the details of the analysis; the results are presented 
and discussed in section 6; section 7 is devoted to the systematic uncertainties; the 
multiplicities measured in this analysis are compared with the results of previous 
experiments in section 8; summary and conclusions are given in section 9. 



2 The experimental setup 

2. 1 The NOMAD detector 

The NOMAD detector (see Fig. 1) has been described in details elsewhere [32]. 
It consists of a number of sub-detectors most of which are located inside a 0.4 T 
dipole magnet with an inner volume of 7.5x3.5x3.5 m^: an active target made of 
drift chambers (DC) [33] with a mass of 2.7 tons, an average density of 0.1 g/cm^ 
and a total length of about one radiation length (~ I.OXq), followed by a tran- 
sition radiation detector (TRD) [34], a preshower detector and an electromagnetic 
calorimeter (ECAL) consisting of 875 lead glass cells [35]. The low density and the 
good instrumentation of the active target (there is less than 1% of a radiation length 
between two consecutive measurements) allow to obtain detailed information about 
the final state of neutrino interactions. The momentum resolution for reconstructed 
tracks is, on average, 3.5% for momenta below 20 GeV/c [32], relevant to most of 
the final state hadrons. 

The TRD is used for electron identification. It provides a pion rejection fac- 
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Fig. 1. Side view of the NOMAD detector. The origin of the coordinate system is on the 
upstream face of the active target. 

tor greater than lO'^ for a 90% electron efficiency in the momentum range from 1 
GeV/c to 50 GeV/c. 

The ECAL provides an energy resolution of 'i.2%/ -/E\G^7] © 1% for elec- 
tromagnetic showers and is essential to measure the total energy flow in neutrino 
interactions. The preshower detector is used for electron and gamma identification. 

A hadron calorimeter is located just after the magnet coil. It is followed by an 
iron absorber and a set of muon chambers. The muon identification efficiency is 
97% for muons of momentum higher than 5 GeV/c. 

The neutrino interaction trigger [36] consists of a coincidence between signals 
from two planes of counters located after the active target in the absence of a signal 
from a large area system of veto counters located upstream of the NOMAD detec- 
tor. 

The NOMAD target consists mainly of carbon (on which about 80% of the 
y^CC interactions occur) and of other elements with similar atomic numbers (mostly 
oxygen and nitrogen). 



2.2 The neutrino beam 



The NOMAD detector is located in the WANF neutrino beam line [37] at 
CERN. In the main focusing mode (positive focusing) the beam is an essentially 
pure i/^ beam (about 93%). In negative focusing it is mainly a z/^ beam of lower pu- 
rity (about 71%). More details about the beam composition can be found in Tables 
1 [38] and 2. 



6 



Neutrino 
flavours 


Flux 


CC interactions in NOMAD 


<E^> [GeV] 


rel.abund. 


<E^> [GeV] 


rel.abund. 




23.5 


1 


43.8 


1 




19.2 


0.0612 


42.8 


0.0255 




37.1 


0.0094 


58.3 


0.0148 




31.3 


0.0024 


54.5 


0.0016 



Table 1 

The CERN SPS neutrino beam composition in the positive focusing mode. 



Neutrino 
flavours 


Flux 


CC interactions in NOMAD 


<E^> [GeV] 


rel.abund. 


<E^> [GeV] 


rel.abund. 




23.5 


0.156 


58 


0.383 




19.6 


1 


33.9 


1 




31.1 


0.0053 


56.3 


0.016 




29.5 


0.0086 


47.8 


0.013 



Table 2 

The CERN SPS neutrino beam composition in the negative focusing mode. 
2.3 The Monte Carlo simulations 



In this analysis the IVIonte Carlo (IVIC) simulation software is based on the 
LEPT06.1 [39] - JETSET7.4 [40] event generator and a full GEANT321 [41] de- 
tector simulation. 

The simulated events passed through the same reconstruction program and se- 
lection criteria that were used for the data. 

An important JETSET parameter for the fragmentation, the fragmentation cut- 
oflFl was set to 0.2 GeV. The strangeness suppression factoiT^ was set to 0.21. 
These values were obtained from a study of hadronic states [42] with invariant 
masses similar to those studied in this analysis. This was also checked by compar- 
ing the fragmentation function of pions and K'g obtained from the IVIC with the one 
obtained from the NOIVIAD data. 

In order to generate /o(980) and /2(1270) with rates comparable to those ob- 
served in our data we changed the default values of the following JETSET parame- 
ters (the default values for both parameters in JETSET7.4 is zero): the probabilit)[l] 
to produce mesons with orbital angular momentum 1 and total angular momentum 



2 PARJ(33) in JETSET7.4 

3 PARJ(2) in JETSET7.4 

4 PARJ(15)inJETSET7.4 
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equal to 0.08; the probabilitjlfj to produce mesons with orbital angular momentum 
1 and total angular momentum 2 equal to 0.185. Note that setting these parameters 
different from zero decreases the number of p°(770) mesons in the MC by the sum 
of values given above (in our case by 26.5%). 

The MC includes a simulation of nuclear effects [43] which takes into account 
the rescattering and absorption of particles inside nuclei. This allowed us to de- 
convolute these effects from our results and extract the resonance yields on free 
nucleons. At our energies and target composition these effects are small. 

The Monte Carlo allows for a direct test of the fitting procedure used to ex- 
tract the number of resonances from the invariant mass plots since the number of 
simulated p°(770) mesons is known. In addition, systematic effects due to other 
resonances can be studied since for most of the reconstructed tracks it is possible 
to match a MC particle at the level of the event generator. 



3 Definition of the kinematical variables 



The following variables are used in this analysis: 
I. The global event variables: 

• Eyis - the total visible energy of the event calculated as the sum of the energies 
of all tracks assigned to the primary vertex and of ECAL clusters identified as 
photons. The reconstructed charged particles were assigned the mass of the pion 
except for identified muons and electrons which were assigned their respective 
masses; 

• AE^is - the measurement error on E^is, calculated from the uncertainty on the 
fitted momentum for each track and from the known energy resolution of the 
ECAL; 

• p™*^* - the modulus of the transverse component of the total event momentum 
with respect to the neutrino beam direction. Lost and mismeasured particles are 
the main contributors to the value of pj^^^^-^ 

• E^ - The neutrino energy in the laboratory system. Monte Carlo studies show 
that the reconstructed hadronic energy is, on the average, underestimated, and 
that this underestimate is correlated with p™*'^^. We have therefore corrected E^s 
by adding a term depending on p^^^'^ and the hadronic energy. E^, is determined 
with a resolution of about 10% and with a systematic error not exceeding 1% in 
the whole energy range of interest; 

• V = Ej, — E^, where E^^^ is the energy of the track identified as a muon: energy 
transfer from the neutrino to the hadrons in the laboratory system. If there were 
more than one muon candidate, the one with the greater transverse momentum 
was chosen; 

5 PARJ(17)inJETSET7.4 
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• ? = P;/ — - the 4-momentum transfer; 

• Q"^ = —q^ - the 4-momentum transfer squared; 

• W'^ = 171% + 2mNV — - the hadronic invariant mass squared {ruN is the 
nucleon mass). 

n. The resonance candidate variables: 

• xp = 2p\\/W {p\\ is the longitudinal momentum of the resonance candidate in 
the rest frame of the hadronic system) - Feynman x variable; 

• z — Elv (where E is the energy of the resonance candidate) - fragmentation 
function variable; 

• Pi_- modulus of the transverse momentum of the resonance candidate relative to 
the total momentum of the hadron jet. 



4 Data sample and event selection 

The present analysis is based on the full NOMAD data sample. For the 
analysis we used the data taken with the neutrino beam set to positive focusing. A 
study of in the positive focusing mode would require additional cuts to reduce 
the admixture of misidentified i/^, which could render a comparison with other 
experiments difficult. For this reason, for the study only the negative focusing 
mode was used. 

A Un{uij)CC event was selected for further processing if 

• there was an identified at the primary vertex with momentum >3 
GeV/c (this cut ensures a good muon identification and a consistency with pre- 
vious neutrino experiments); 

• the invariant mass of the hadronic system W was higher than 2 GeV (in order to 
limit the study to the deep inelastic region); 

• there were at least 3 tracks (including the muon) at the primary vertex (this is 
necessary to build a resonance candidate). 

The following quality cuts were applied in order to ensure an accurate event 
reconstruction: 

• the primary vertex (see Fig. 1 for the coordinate system) should be inside the 
following fiducial volume: 

\X, Y\ < 120cm, 10cm< Z < 390cm; 

• the event should be measured with reasonable accuracy (AE^is/E^ig < 30%). 

After cuts we obtain 669252 v,.CC and 15927 V,.CC event candidates. 
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5 Extraction of the resonance signal 



All combinations of tracks with momenta larger than 0.1 GeV/c originating 
at the primary vertex and not identified as electrons or muons were used for the 
construction of the resonance candidates. All used tracks were assigned the pion 
mass. 
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Fig. 2. Raw distribution of the tt'^it in- 
variant mass in the v^CC sample. 



Fig. 3. Raw distribution of the invariant 
mass of TT'^TT^" and 7r~7r~ pairs in the 
i^i_iCC sample. The line is the result of 
the fit described in the text. 



o 
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Fig. 4. Raw distribution of the tt'^tt invariant mass with xp > 0.6. 
The distribution of the 71~^ti~ invariant mass is shown in Fig. 2. For comparison 
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Fig. 5. TT+vr invariant mass distribution of Monte Carlo events including the contributions 
of signal, reflections and combinatorial background. 

the distribution of the invariant mass of tt+tt"^ and 7r~7r~ pairs is shown in Fig. 3. 
The TT+TT^ distribution shows an enhancement at the p'^(770) mass which is not 
present in the like-sign distribution. When a cut xp > 0.6 is applied to reduce the 
combinatorial background, clear peaks at the /o(980) and /2(1270) masses become 
visible, as shown in Fig. 4. 

We now address the treatment of reflections. Reflections are contributions to 
the TT+TT" invariant mass distribution from resonances decaying to tt+tt" + some 
other particles (e.g. cu — > n^n^n^) or contributions from resonances decaying to 
particles other than pions and thus assigned wrong mass values. The reflections 
were estimated from the fully simulated MC events and taken into account in the fit 
of the TT+TT" invariant mass distribution (Fig. 5). It is possible however that the mul- 
tiplicity of the resonances giving reflections in the data could be different from the 
one in the MC. In order to take this into account the sum of all reflections obtained 
from the MC was multiplied by an overall normalization factor p which was left 
as a free parameter in the fit. In NOMAD an additional contribution to the tt+tt" 
invariant mass distribution arises from early decays of Kg and misidentified 7 con- 
versions (see Fig. 5). For distances less than 15cm from the primary vertex along 
the beam axis it is difficult to distinguish tracks beginning at the primary vertex 
from tracks originating at a vertex. This contribution was also estimated from 
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the MC simulation and taken into account together with other reflections. As can 
be seen in Fig. 5, the largest contribution is from uj vt+tt'tt^. 

Ignoring reflections could result in incorrect determination of the resonance pa- 
rameters. As an example, in our case this would result in a p° mass value extracted 
from the fit about 10 MeV lower than the accepted value [44]. 

The resonance signal is determined by fitting the invariant mass distribution of 
all possible tt+tt" combinations, dN/ dm, to the expression 

dN 

— ^[1 + BW{m)]BG{m)+p- RF{m), (1) 



where 

BW{m) = aiBWp{m) + a2BWf^{m) + asBWf^im) (2) 



is a relativistic Breit-Wigner function [45], BG(m) is the combinatorial background 
and RF{m) is the contribution from reflections. 
The Breit-Wigner function is 



BW = - (3) 



where 



^R-^nih''^'— (4) 
kn m 



Here mu and F^ are the mass and width of the resonance R; L is the relative 
orbital angular momentum of the two pions (equal to the spin of the resonance): 
L=0 for /o(980), L=l for p°(770), and L=2 for /2(1270); k is the pion momentum 
in the resonance rest frame; is the value of k when m = uir. 

The background was assumed to have the following shape 

BG = 04(777. — 2m,T^)"'^ exp{aQ'm + a'jvr? + ogm^) (5) 



which takes into account the threshold effect and the exponential fall-off of the dis- 
tribution at high values of m. 

In Equation (1) the Breit-Wigner functions are multiplied by BG in order to 
account for the available phase-space. The parameters a\ to as, v and the masses 
and widths of the three resonances are free fit parameters. The mass range covered 
by the fit was 0.3 GeV< m <2 GeV. 

As a cross-check of the reliability of the procedure we performed the fit using 
another parameterization suggested by A.Minaenko [46]. The p°(770) yield found 
by this method agrees with our result within 0.5%. 



12 



A fit to the distribution of the invariant mass for Uke sign pairs (Fig. 3) gives 
a resonance yield at niR — 770MeV compatible with zero within one standard 
deviation. 

In order to take into account the experimental resolution on the tt+tt" invariant 
mass we followed the method described in [9]. We estimated the invariant mass 
resolution function from the errors on the track parameters in the data as given by 
the reconstruction program. This function has a FWHM of about 25 MeV and it 
was convoluted with expression (1) for the fit. 

The procedure used to obtain the results presented in section 6 proceeded 
through the steps listed below. Steps (i) and (ii) refer to the procedure of obtain- 
ing the overall average multiplicities and the resonance widths and masses, while 
the subsequent steps were made to obtain the differential distributions for the kine- 
matical variables. 

(1) Expression (1) is fitted to the raw inclusive experimental distribution. 

(2) The "fully simulated" MC is used to estimate the reconstruction efficiency of 
various resonances by comparing the yields of reconstructed and generated 
resonances. The uncorrected yields measured in the data are then divided by 
these efficiencies. 

(3) A fit to the experimental distributions in each bin of kinematical variables is 
performed. The reflection normalization factor p and the masses and widths 
of the resonances are fixed at the values found in step (i). 

(4) Correction factors are obtained in each of the kinematical variable bins sim- 
ilarly to step (ii). The yields of the resonances in each bin are corrected ac- 
cordingly. 

(5) The corrected experimental distributions are plotted and compared with the 
MC distributions at the event generator level. 

The P^CC sample was treated in the same way, however, due to its lower statis- 
tics we limited the analysis to the steps (i) - (ii) and to the p°(770) resonance only. 



6 The analysis 

6. 1 Resonance multiplicities and properties 

All the results reported in this section are obtained from the analysis of the 
v^CC sample. 

In Fig. 6 we show the data after subtracting the contribution of the reflections 
as obtained from the fit. In Fig. 7 we further subtract the contribution of the com- 
binatorial background. The three resonances are clearly visible. 

After correcting the result of the fit for the efficiencies as determined from 
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the MC (about 0.77 for all the studied resonances) and for the branching ratios 
BR{R — > 7r+7r~) [44] we obtain the results shown in Table 3. 



Resonance 


Branching 
ratio TT+vr^ [44] 


Number of 
resonances 


Average 
Multiplicity 


Mass 
(MeV) 


r(MeV) 


p0(770) 


1.000 


130368±4336 


0.195±0.007 


768±2 


151±7 


/o(980) 


0.666 


11809±1965 


0.018±0.003 


963±5 


35±10 


/2(1270) 


0.564 


25189±3958 


0.038±0.006 


1286±9 


198±30 



Table 3 



Corrected numbers and multiplicities of the three resonances and their measured masses 
and widths. The errors are statistical only. 




Fig. 6. Distribution of the tt'^tt" invari- 
ant mass after subtracting the contribu- 
tion of reflections. 



Fig. 7. Distribution of the tt'^tt" invari- 
ant mass after subtracting the contribu- 
tions of both reflections and combinato- 
rial background. 
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6.2 Resonance multiplicities as a function of the kinematical variables 



These results were obtained as described in the section 5, steps (iii) - (v). For 
all bins of all kinematical distributions the resonance signals were clearly visible 
and the value of the reduced was around 1.0. All the numbers corrected for 
experimental efficiencies and branching ratios, are compared with the results of the 
modified Lund simulation (see section 2.3). The errors shown in the figures are 
statistical only. 

The increase of the average multiplicity of the three resonances as a function 




10 10^ 10 10^ 10 10^ 

E,(GeV) E.(GeV) E,(GeV) 

Fig. 8. Average p°(770)('aj, fo{980)(b) and /2(1270)('cJ multiplicity as a function of Ey. 
The solid line represents the result of the modified Lund simulation. The errors are statisti- 
cal only. 

of E^, is shown in Fig. 8. 




Fig. 9. Average p°(770)('aj, fQ{dm)(b) and /2(1270)('cj multiplicity as a function of W . 
The solid line represents the result of the modified Lund simulation. The dashed line repre- 
sents a fitted function a + 6 x InW^. The errors are statistical only. 

Fig. 9 shows the average multiplicities as a function of W"^. A fit to the average 
multiplicities with the expression a + 6 x InW"^ gives the result presented in Table 4 
and shown as a dashed line in Fig. 9. 
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Resonance 


a 


b 


p0(770) 


0.026 ± 0.024 


0.062 ± 0.009 


/o(980) 


-0.018 ± 0.008 


0.038 ± 0.009 


/2(1270) 


-0.011 ± 0.003 


0.021 ± 0.004 



Table 4 

The parameters of the W dependence of the average multiplicities fitted to the expression 
a + 6 X InW'^. The errors are statistical only. 

The dependence of the average multiplicities on is plotted in Fig. 10. 




1 10 1 10 

(GeV^) (GeV^) (GeV^) 



Fig. 10. Average p^(770)(a), fo(980)(b) and /2(1270)('cj multiplicity as a function ofQ"^. 
The solid line represents the result of the Lund simulation. 

Fig. 11 shows the xp dependence of the resonance production. From this pic- 
ture we see that the forward production of p°(770) and /2(1270) is enhanced with 
respect to the backward production. The Lund model overestimates the backward 
production of these resonances. 

The resonance fragmentation function is shown in Fig. 12. The high z region 




(b)W980) 



(c)f,(1270) 




-1 -0.8-0.6-0.4-0.2 0.2 0.4 0.6 0.8 1 -1 -0.8-0.6-0.4-0.2 0.2 0.4 0.6 0.8 1 -1 -0.8-0.6-0.4-0.2 0.2 0.4 0.6 0.8 1 



Fig. 11. Average p°(770)('flj, /o(980)ffoj and /2(1270)('cj multiplicity as a function ofxp- 
The solid line represents the result of the Lund simulation 

of the distributions is overestimated by the Lund model. 

The p\ dependence of the resonance production is shown in Fig. 13. In the 
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case of p°(770) our result favours a steeper p\ slope than previously measured in 
[18,20] and agrees with that of ref. [23]. A fit using the function ae~^^ gives 

a = (1.28 ± 0.04)(GeV/c)-2, b = (4.79 ± 0.07)(GeV/c)-2. 

An experiment [18] found that the slope parameter 6 in the range > 0.5(GeV/c)^ 
is different from the slope below this value. We find 

6 = (5.3±0.2)(GeV/c)-2. 

for pi < 0.5(GeV/c)2; and 

6 = (4.1±0.1)(GeV/c)-2. 

for p']_ > 0.5(GeV/c)^. Thus we confirm the existence of two different slopes as 
found earlier [18]. 

We also observe that the distribution for /2(1270) in the data is harder than 
the results of the Lund model simulation. This could suggest a mass dependence of 




Fig. 12. Average p"(770)faj, fo{980)(b) and f2{l270)(c) multiplicity as a function of z 
(fragmentation functions). The solid line represents the result of the Lund simulation. 




Fig. \3. Average p^{770)(a), fo{980)(b), and f2{1270)(c) multiplicity as a function of p,. 
The solid line represents the result of the Lund simulation. The dashed straight line repre- 
sents the result of the fit by the function ae~^-i-. 
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the p\ behaviour. 

Finally, we address the possibility that the /o(980)(b) corresponds to the vac- 
uum scalar state suggested by Gribov (see Section 1). Fig. 14 shows the resonance 
multiplicities divided by the ones simulated with the Lund model as a function 
of the hadron jet charged track multiplicity. The measured and MC simulated reso- 
nance multiplicities were normalized to the same average multiplicity. From Fig. 14 
we conclude that in our W range and at our level of accuracy (15 - 20%) we do not 
observe the /o(980) excess at low multiplicities predicted by this model. 




2 4 6 8 10 12 2 4 6 8 10 12 2 4 6 8 10 12 

"ch "ch "ch 

Fig. 14. The resonance multiplicities, normalized to the ones simulated with the Lund model 
as a function of the hadron jet charged particle multiplicity. 



7 Systematic uncertainties 

The main source of systematic uncertainty results from the subtraction of the 
reflections. Reflections in the data can be different from the MC not only by a 
single overall factor but could also have a different shape due to different relative 
resonance yields. We estimated this effect as follows. The fit to the vr+vr" invari- 
ant mass distribution gave the reflection normalization factor p. We repeated the 
fit changing p beyond the statistical uncertainty obtained from the initial fit (up to 
15%). This results in relative changes of the resonances yields as reported in Ta- 
ble 5. 

Other sources of systematics are: 

• The uncertainty due to the choice of the fit boundaries: increasing the left bound- 
ary by 25 MeV or lowering the right boundary by 300 MeV changes the yields 

as shown in Table 5; 

• Removing the cubic or square term in the exponential function used to describe 
the background changes the result as shown in the line labelled "Fit model" of 
Table 5; 

• Changing the bin width has a negligible effect on the result. 
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Source of systematics 


p0(770) 


/o(980) 


/2(1270) 


Reflection subtraction 


8% 


10% 


9% 


Fit boundaries 


4% 


8% 


9% 


Fit model 


2% 


9% 


13% 


Total 


9% 


16% 


18% 



Table 5 

Systematic uncertainties on the resonance yields. 

The overall systematic uncertainty on the measured resonance yields, also shown 
in Table 5, is obtained by adding all effects in quadrature. 

The systematic uncertainties affecting the resonance widths are shown in Ta- 
ble 6. 



Source of systematics 


/(770) 


/o(980) 


/2(1270) 


Reflection subtraction 


13 


3 


11 


Bin width 


5 


3 


1 


Fit model 


12 


5 


6 


Total 


18 


7 


13 



Table 6 

Systematic uncertainties on the resonance widths, in MeV. 

Table 7 summarizes our results with statistical and systematic uncertainties 
added in quadrature. 



Resonance 


Number of 

resonances 


Average 
Multiplicity 


Mass 
(MeV) 


r(MeV) 


p0(770) 


130368±12509 


0.195±0.019 


768±2.5 


151±19 


/o(980) 


11809±2726 


0.018±0.004 


963±5 


35±12 


/2(1270) 


25189±6018 


0.038±0.009 


1286±11 


198±32 



Table 7 

Total corrected numbers and yields of resonances and their masses and widths. The errors 
include statistical and systematic error. 

8 Yields of p°(770) and comparison with previous experiments 

The average p°(770) multiplicity measured in the NOMAD experiment for neutrino 

and antineutrino charged current interactions is shown in Table 8, together with 
the results from previous experiments. In the neutrino charged current sample, our 
measurement confirms the results of all previous experiments in a similar W range 
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with a much better statistical accuracy. The deviation of the SKAT [24] result is 
usually explained as due to the lower average W. In the case of z/^, there was 
some disagreement between previous experiments which could not be explained 
by different W ranges. Our result confirms the measurements in which the average 
p°(770) multiplicity in V^CC interactions is the same as in v^CC interactions. 



Reaction 


W 


Average multiplicity 


Reference 


f^nucleon 


W>2GeV 


0.195±0.007±0.019 


this analysis 




5.5 


0.18±0.02±0.02 


[18] 


z^mP 


W>2GeV 


0.21±0.03±0.03 


[23] 




4.9 


0.16±0.02 


[21] 


z^mP 


W>2GeV 


0.20±0.05 


[28] 


z^mP 


5.0 


0.21 ±0.04 


[25] 


z^^n 


W>2GeV 


0.21±0.02±0.03 


[23] 


z^m'i 




0.17±0.05 


[28] 




4.8 


0.17±0.04 


[19] 


fjuFreon 


3.1 


0.09±0.02 


[24] 


I7^nucleon 


W>2GeV 


0.18±0.02±0.02 


this analysis 


^nV 


4.4 


0.12±0.02±0.02 


[18] 


z^mP 


W>2GeV 


0.20±0.02±0.03 


[23] 


z^mP 


4.2 


0.11 ±0.02 


[22] 


Z^AtP 


3.4 


0.21 ±0.03 


[26] 


Z^luH 


W>2GeV 


0.18±0.03±0.03 


[23] 




3.9 


0.12±0.02 


[20] 


I7^Ne 


4.3 


0.13±0.02 


[27] 



Table 8 

Average p^{770) multiplicity in andV^ experiments. 
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9 Conclusions 



The inclusive production of the meson resonance p°(770) in neutrino - nucleon 
charged current interactions has been studied with much better statistical accuracy 
than in previous experiments. 

We also looked for other resonances decaying to tt+tt" pairs. For the first time 
the /o(980) meson is observed in neutrino interactions. The signal has a signifi- 
cance of about 5 standard deviations. 

We see a clear signal of /2( 12 70). Its production in neutrino interactions is reli- 
ably established. 

A global fit to the tt+tt" pair mass distribution allows us to obtain the masses 
and the widths of the p°(770), /o(980) and /2(1270). 

The average multiplicities of p°(770), /o(980) and /2(1270) are measured as a 
function of VF^, and other kinematic variables. These results are compared with 
a simulation based on a modified Lund model. Good agreement for most of the 
distributions is found. 

In addition, the production of p°(770) in CC interactions has been studied 
with lower statistics. This allowed us to clarify the ambiguous experimental situa- 
tion on the production of p°(770) in antineutrino interactions. 

We see no evidence of the enhanced /o(980) production in low multiplicity 
events predicted by Gribov et al. [14], [15]. 

These studies demonstrate the capability of NOMAD to obtain results which up 
to now have been the exclusive domain of bubble-chamber experiments, but with 
much improved statistics. 
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